Gene-educational attainment interactions in a multi-population genome-wide meta-analysis identify novel lipid loci.

Introduction: Educational attainment, widely used in epidemiologic studies as a surrogate for socioeconomic status, is a predictor of cardiovascular health outcomes. Methods: A two-stage genome-wide meta-analysis of low-density lipoprotein cholesterol (LDL), high-density lipoprotein cholesterol (HDL), and triglyceride (TG) levels was performed while accounting for gene-educational attainment interactions in up to 226,315 individuals from five population groups. We considered two educational attainment variables: "Some College" (yes/no, for any education beyond high school) and "Graduated College" (yes/no, for completing a 4-year college degree). Genome-wide significant (p < 5 × 10-8) and suggestive (p < 1 × 10-6) variants were identified in Stage 1 (in up to 108,784 individuals) through genome-wide analysis, and those variants were followed up in Stage 2 studies (in up to 117,531 individuals). Results: In combined analysis of Stages 1 and 2, we identified 18 novel lipid loci (nine for LDL, seven for HDL, and two for TG) by two degree-of-freedom (2 DF) joint tests of main and interaction effects. Four loci showed significant interaction with educational attainment. Two loci were significant only in cross-population analyses. Several loci include genes with known or suggested roles in adipose (FOXP1, MBOAT4, SKP2, STIM1, STX4), brain (BRI3, FILIP1, FOXP1, LINC00290, LMTK2, MBOAT4, MYO6, SENP6, SRGAP3, STIM1, TMEM167A, TMEM30A), and liver (BRI3, FOXP1) biology, highlighting the potential importance of brain-adipose-liver communication in the regulation of lipid metabolism. An investigation of the potential druggability of genes in identified loci resulted in five gene targets shown to interact with drugs approved by the Food and Drug Administration, including genes with roles in adipose and brain tissue. Discussion: Genome-wide interaction analysis of educational attainment identified novel lipid loci not previously detected by analyses limited to main genetic effects.

Atherosclerotic cardiovascular disease risk and small dense low-density lipoprotein cholesterol in men, women, African Americans and non-African Americans: The pooling project.

BACKGROUND AND AIMS: Elevated small dense low-density lipoprotein-cholesterol (sdLDL-C) has been reported to be associated with increased atherosclerotic cardiovascular disease (ASCVD) risk. Our aims were to determine whether direct and calculated sdLDL-C were significant independent ASCVD risk factors in sex and race subgroups. METHODS: In a total of 15,933 participants free of ASCVD at baseline (median age 62 years, 56.7% female, 19.7% African American) fasting plasma lipids and sdLDL-C were measured by standardized automated methods. All subjects were followed for 10 years for incident ASCVD, which developed in 9.7% of subjects. SdLDL-C values were also calculated. Univariate and multivariate analyses were carried out to assess for independent associations with incident ASCVD after adjustment for all standard risk factors. RESULTS: All standard risk factors were significantly associated with incident ASCVD on univariate analysis, as were direct and calculated sdLDL-C. These latter parameters were also significant when added to the pooled cohort risk equation. However, associations were significantly stronger for direct sdLDL-C and were not significant for calculated values once direct values were in the model. In contrast to calculated values, top quartile direct sdLDL-C was significantly independently associated with incident ASCVD versus bottom quartile values in all subjects and subgroups, except African Americans (hazard ratios ≥1.50, p < 0.01). Subjects with direct values ≥ 50 mg/dL versus <25 mg/dL had significantly higher independent ASCVD risk in all groups (hazard ratios >1.49, all p < 0.01). CONCLUSIONS: Having a direct small dense low-density lipoprotein cholesterol value ≥ 50 mg/dL is a significant independent ASCVD risk-enhancer.

Whole-exome sequence analysis of anthropometric traits illustrates challenges in identifying effects of rare genetic variants.

Anthropometric traits, measuring body size and shape, are highly heritable and significant clinical risk factors for cardiometabolic disorders. These traits have been extensively studied in genome-wide association studies (GWASs), with hundreds of genome-wide significant loci identified. We performed a whole-exome sequence analysis of the genetics of height, body mass index (BMI) and waist/hip ratio (WHR). We meta-analyzed single-variant and gene-based associations of whole-exome sequence variation with height, BMI, and WHR in up to 22,004 individuals, and we assessed replication of our findings in up to 16,418 individuals from 10 independent cohorts from Trans-Omics for Precision Medicine (TOPMed). We identified four trait associations with single-nucleotide variants (SNVs; two for height and two for BMI) and replicated the LECT2 gene association with height. Our expression quantitative trait locus (eQTL) analysis within previously reported GWAS loci implicated CEP63 and RFT1 as potential functional genes for known height loci. We further assessed enrichment of SNVs, which were monogenic or syndromic variants within loci associated with our three traits. This led to the significant enrichment results for height, whereas we observed no Bonferroni-corrected significance for all SNVs. With a sample size of ∼20,000 whole-exome sequences in our discovery dataset, our findings demonstrate the importance of genomic sequencing in genetic association studies, yet they also illustrate the challenges in identifying effects of rare genetic variants.

Exome sequence association study of levels and longitudinal change of cardiovascular risk factor phenotypes in European Americans and African Americans from the Atherosclerosis Risk in Communities Study.

Cardiovascular disease (CVD) is responsible for 31% of all deaths worldwide. Among CVD risk factors are age, race, increased systolic blood pressure (BP), and dyslipidemia. Both BP and blood lipids levels change with age, with a dose-dependent relationship between the cumulative exposure to hyperlipidemia and the risk of CVD. We performed an exome sequence association study using longitudinal data with up to 7805 European Americans (EAs) and 3171 African Americans (AAs) from the Atherosclerosis Risk in Communities (ARIC) study. We assessed associations of common (minor allele frequency > 5%) nonsynonymous and splice-site variants and gene-based sets of rare variants with levels and with longitudinal change of seven CVD risk factor phenotypes (BP traits: systolic BP, diastolic BP, pulse pressure; lipids traits: triglycerides, total cholesterol, high-density lipoprotein cholesterol [HDL-C], low-density lipoprotein cholesterol [LDL-C]). Furthermore, we investigated the relationship of the identified variants and genes with select CVD endpoints. We identified two novel genes: DCLK3 associated with the change of HDL-C levels in AAs and RAB7L1 associated with the change of LDL-C levels in EAs. RAB7L1 is further associated with an increased risk of heart failure in ARIC EAs. Investigation of the contribution of genetic factors to the longitudinal change of CVD risk factor phenotypes promotes our understanding of the etiology of CVD outcomes, stressing the importance of incorporating the longitudinal structure of the cohort data in future analyses.

Multi-ancestry genome-wide gene-sleep interactions identify novel loci for blood pressure.

Long and short sleep duration are associated with elevated blood pressure (BP), possibly through effects on molecular pathways that influence neuroendocrine and vascular systems. To gain new insights into the genetic basis of sleep-related BP variation, we performed genome-wide gene by short or long sleep duration interaction analyses on four BP traits (systolic BP, diastolic BP, mean arterial pressure, and pulse pressure) across five ancestry groups in two stages using 2 degree of freedom (df) joint test followed by 1df test of interaction effects. Primary multi-ancestry analysis in 62,969 individuals in stage 1 identified three novel gene by sleep interactions that were replicated in an additional 59,296 individuals in stage 2 (stage 1 + 2 Pjoint < 5 × 10-8), including rs7955964 (FIGNL2/ANKRD33) that increases BP among long sleepers, and rs73493041 (SNORA26/C9orf170) and rs10406644 (KCTD15/LSM14A) that increase BP among short sleepers (Pint < 5 × 10-8). Secondary ancestry-specific analysis identified another novel gene by long sleep interaction at rs111887471 (TRPC3/KIAA1109) in individuals of African ancestry (Pint = 2 × 10-6). Combined stage 1 and 2 analyses additionally identified significant gene by long sleep interactions at 10 loci including MKLN1 and RGL3/ELAVL3 previously associated with BP, and significant gene by short sleep interactions at 10 loci including C2orf43 previously associated with BP (Pint < 10-3). 2df test also identified novel loci for BP after modeling sleep that has known functions in sleep-wake regulation, nervous and cardiometabolic systems. This study indicates that sleep and primary mechanisms regulating BP may interact to elevate BP level, suggesting novel insights into sleep-related BP regulation.

Discovery and fine-mapping of height loci via high-density imputation of GWASs in individuals of African ancestry.

Although many loci have been associated with height in European ancestry populations, very few have been identified in African ancestry individuals. Furthermore, many of the known loci have yet to be generalized to and fine-mapped within a large-scale African ancestry sample. We performed sex-combined and sex-stratified meta-analyses in up to 52,764 individuals with height and genome-wide genotyping data from the African Ancestry Anthropometry Genetics Consortium (AAAGC). We additionally combined our African ancestry meta-analysis results with published European genome-wide association study (GWAS) data. In the African ancestry analyses, we identified three novel loci (SLC4A3, NCOA2, ECD/FAM149B1) in sex-combined results and two loci (CRB1, KLF6) in women only. In the African plus European sex-combined GWAS, we identified an additional three novel loci (RCCD1, G6PC3, CEP95) which were equally driven by AAAGC and European results. Among 39 genome-wide significant signals at known loci, conditioning index SNPs from European studies identified 20 secondary signals. Two of the 20 new secondary signals and none of the 8 novel loci had minor allele frequencies (MAF) < 5%. Of 802 known European height signals, 643 displayed directionally consistent associations with height, of which 205 were nominally significant (p < 0.05) in the African ancestry sex-combined sample. Furthermore, 148 of 241 loci contained ≤20 variants in the credible sets that jointly account for 99% of the posterior probability of driving the associations. In summary, trans-ethnic meta-analyses revealed novel signals and further improved fine-mapping of putative causal variants in loci shared between African and European ancestry populations.

Small Dense Low-Density Lipoprotein Cholesterol Is the Most Atherogenic Lipoprotein Parameter in the Prospective Framingham Offspring Study.

Background Elevated plasma levels of direct low-density lipoprotein cholesterol (LDL-C), small dense LDL-C (sdLDL-C), low-density lipoprotein (LDL) triglycerides, triglycerides, triglyceride-rich lipoprotein cholesterol, remnant lipoprotein particle cholesterol, and lipoprotein(a) have all been associated with incident atherosclerotic cardiovascular disease (ASCVD). Our goal was to assess which parameters were most strongly associated with ASCVD risk. Methods and Results Plasma total cholesterol, triglycerides, high-density lipoprotein cholesterol, direct LDL-C, sdLDL-C, LDL triglycerides, remnant lipoprotein particle cholesterol, triglyceride-rich lipoprotein cholesterol, and lipoprotein(a) were measured using standardized automated analysis (coefficients of variation, <5.0%) in samples from 3094 fasting subjects free of ASCVD. Of these subjects, 20.2% developed ASCVD over 16 years. On univariate analysis, all ASCVD risk factors were significantly associated with incident ASCVD, as well as the following specialized lipoprotein parameters: sdLDL-C, LDL triglycerides, triglycerides, triglyceride-rich lipoprotein cholesterol, remnant lipoprotein particle cholesterol, and direct LDL-C. Only sdLDL-C, direct LDL-C, and lipoprotein(a) were significant on multivariate analysis and net reclassification after adjustment for standard risk factors (age, sex, hypertension, diabetes mellitus, smoking, total cholesterol, and high-density lipoprotein cholesterol). Using the pooled cohort equation, many specialized lipoprotein parameters individually added significant information, but no parameter added significant information once sdLDL-C (hazard ratio, 1.42; P<0.0001) was in the model. These results for sdLDL-C were confirmed by adjusted discordance analysis versus calculated non-high-density lipoprotein cholesterol, in contrast to LDL triglycerides. Conclusions sdLDL-C, direct LDL-C, and lipoprotein(a) all contributed significantly to ASCVD risk on multivariate analysis, but no parameter added significant risk information to the pooled cohort equation once sdLDL-C was in the model. Our data indicate that small dense LDL is the most atherogenic lipoprotein parameter.

Methylome-wide association study of central adiposity implicates genes involved in immune and endocrine systems.

Aim: We conducted a methylome-wide association study to examine associations between DNA methylation in whole blood and central adiposity and body fat distribution, measured as waist circumference, waist-to-hip ratio and waist-to-height ratio adjusted for body mass index, in 2684 African-American adults in the Atherosclerosis Risk in Communities study. Materials & methods: We validated significantly associated cytosine-phosphate-guanine methylation sites (CpGs) among adults using the Women's Health Initiative and Framingham Heart Study participants (combined n = 5743) and generalized associations in adolescents from The Raine Study (n = 820). Results & conclusion: We identified 11 CpGs that were robustly associated with one or more central adiposity trait in adults and two in adolescents, including CpG site associations near TXNIP, ADCY7, SREBF1 and RAP1GAP2 that had not previously been associated with obesity-related traits.

Efficient gene-environment interaction tests for large biobank-scale sequencing studies.

Complex human diseases are affected by genetic and environmental risk factors and their interactions. Gene-environment interaction (GEI) tests for aggregate genetic variant sets have been developed in recent years. However, existing statistical methods become rate limiting for large biobank-scale sequencing studies with correlated samples. We propose efficient Mixed-model Association tests for GEne-Environment interactions (MAGEE), for testing GEI between an aggregate variant set and environmental exposures on quantitative and binary traits in large-scale sequencing studies with related individuals. Joint tests for the aggregate genetic main effects and GEI effects are also developed. A null generalized linear mixed model adjusting for covariates but without any genetic effects is fit only once in a whole genome GEI analysis, thereby vastly reducing the overall computational burden. Score tests for variant sets are performed as a combination of genetic burden and variance component tests by accounting for the genetic main effects using matrix projections. The computational complexity is dramatically reduced in a whole genome GEI analysis, which makes MAGEE scalable to hundreds of thousands of individuals. We applied MAGEE to the exome sequencing data of 41,144 related individuals from the UK Biobank, and the analysis of 18,970 protein coding genes finished within 10.4 CPU hours.

Role of Rare and Low-Frequency Variants in Gene-Alcohol Interactions on Plasma Lipid Levels.

BACKGROUND: Alcohol intake influences plasma lipid levels, and such effects may be moderated by genetic variants. We aimed to characterize the role of aggregated rare and low-frequency protein-coding variants in gene by alcohol consumption interactions associated with fasting plasma lipid levels. METHODS: In the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium, fasting plasma triglycerides and high- and low-density lipoprotein cholesterol were measured in 34 153 individuals with European ancestry from 5 discovery studies and 32 277 individuals from 6 replication studies. Rare and low-frequency functional protein-coding variants (minor allele frequency, ≤5%) measured by an exome array were aggregated by genes and evaluated by a gene-environment interaction test and a joint test of genetic main and gene-environment interaction effects. Two dichotomous self-reported alcohol consumption variables, current drinker, defined as any recurrent drinking behavior, and regular drinker, defined as the subset of current drinkers who consume at least 2 drinks per week, were considered. RESULTS: We discovered and replicated 21 gene-lipid associations at 13 known lipid loci through the joint test. Eight loci (PCSK9, LPA, LPL, LIPG, ANGPTL4, APOB, APOC3, and CD300LG) remained significant after conditioning on the common index single-nucleotide polymorphism identified by previous genome-wide association studies, suggesting an independent role for rare and low-frequency variants at these loci. One significant gene-alcohol interaction on triglycerides in a novel locus was significantly discovered (P=6.65×10-6 for the interaction test) and replicated at nominal significance level (P=0.013) in SMC5. CONCLUSIONS: In conclusion, this study applied new gene-based statistical approaches and suggested that rare and low-frequency genetic variants interacted with alcohol consumption on lipid levels.

A unified method for rare variant analysis of gene-environment interactions.

Advanced technology in whole-genome sequencing has offered the opportunity to comprehensively investigate the genetic contribution, particularly rare variants, to complex traits. Several region-based tests have been developed to jointly model the marginal effect of rare variants, but methods to detect gene-environment (GE) interactions are underdeveloped. Identifying the modification effects of environmental factors on genetic risk poses a considerable challenge. To tackle this challenge, we develop a method to detect GE interactions for rare variants using generalized linear mixed effect model. The proposed method can accommodate either binary or continuous traits in related or unrelated samples. Under this model, genetic main effects, GE interactions, and sample relatedness are modeled as random effects. We adopt a kernel-based method to leverage the joint information across rare variants and implement variance component score tests to reduce the computational burden. Our simulation studies of continuous and binary traits show that the proposed method maintains correct type I error rates and appropriate power under various scenarios, such as genotype main effects and GE interaction effects in opposite directions and varying the proportion of causal variants in the model. We apply our method in the Framingham Heart Study to test GE interaction of smoking on body mass index or overweight status and replicate the Cholinergic Receptor Nicotinic Beta 4 gene association reported in previous large consortium meta-analysis of single nucleotide polymorphism-smoking interaction. Our proposed set-based GE test is computationally efficient and is applicable to both binary and continuous phenotypes, while appropriately accounting for familial or cryptic relatedness.

Multi-ancestry sleep-by-SNP interaction analysis in 126,926 individuals reveals lipid loci stratified by sleep duration.

Both short and long sleep are associated with an adverse lipid profile, likely through different biological pathways. To elucidate the biology of sleep-associated adverse lipid profile, we conduct multi-ancestry genome-wide sleep-SNP interaction analyses on three lipid traits (HDL-c, LDL-c and triglycerides). In the total study sample (discovery + replication) of 126,926 individuals from 5 different ancestry groups, when considering either long or short total sleep time interactions in joint analyses, we identify 49 previously unreported lipid loci, and 10 additional previously unreported lipid loci in a restricted sample of European-ancestry cohorts. In addition, we identify new gene-sleep interactions for known lipid loci such as LPL and PCSK9. The previously unreported lipid loci have a modest explained variance in lipid levels: most notable, gene-short-sleep interactions explain 4.25% of the variance in triglyceride level. Collectively, these findings contribute to our understanding of the biological mechanisms involved in sleep-associated adverse lipid profiles.

Direct Versus Calculated LDL Cholesterol and C-Reactive Protein in Cardiovascular Disease Risk Assessment in the Framingham Offspring Study.

BACKGROUND: Increases in circulating LDL cholesterol (LDL-C) and high-sensitivity C-reactive protein (hsCRP) concentrations are significant risk factors for cardiovascular disease (CVD). We assessed direct LDL-C and hsCRP concentrations compared to standard risk factors in the Framingham Offspring Study. METHODS: We used stored frozen plasma samples (-80 °C) obtained after an overnight fast from 3147 male and female participants (mean age, 58 years) free of CVD at cycle 6 of the Framingham Offspring Study. Overall, 677 participants (21.5%) had a CVD end point over a median of 16.0 years of follow-up. Total cholesterol (TC), triglyceride (TG), HDL cholesterol (HDL-C), direct LDL-C (Denka Seiken and Kyowa Medex methods), and hsCRP (Dade Behring method) concentrations were measured by automated analysis. LDL-C was also calculated by both the Friedewald and Martin methods. RESULTS: Considering all CVD outcomes on univariate analysis, significant factors included standard risk factors (age, hypertension, HDL-C, hypertension treatment, sex, diabetes, smoking, and TC concentration) and nonstandard risk factors (non-HDL-C, direct LDL-C and calculated LDL-C, TG, and hsCRP concentrations). On multivariate analysis, only the Denka Seiken direct LDL-C and the Dade Behring hsCRP were still significant on Cox regression analysis and improved the net risk reclassification index, but with modest effects. Discordance analysis confirmed the benefit of the Denka Seiken direct LDL-C method for prospective hard CVD endpoints (new-onset myocardial infarction, stroke, and/or CVD death). CONCLUSIONS: Our data indicate that the Denka Seiken direct LDL-C and Dade Behring hsCRP measurements add significant, but modest, information about CVD risk, compared to standard risk factors and/or calculated LDL-C.

Multi-ancestry genome-wide gene-smoking interaction study of 387,272 individuals identifies new loci associated with serum lipids.

The concentrations of high- and low-density-lipoprotein cholesterol and triglycerides are influenced by smoking, but it is unknown whether genetic associations with lipids may be modified by smoking. We conducted a multi-ancestry genome-wide gene-smoking interaction study in 133,805 individuals with follow-up in an additional 253,467 individuals. Combined meta-analyses identified 13 new loci associated with lipids, some of which were detected only because association differed by smoking status. Additionally, we demonstrate the importance of including diverse populations, particularly in studies of interactions with lifestyle factors, where genomic and lifestyle differences by ancestry may contribute to novel findings.

Network analysis of drug effect on triglyceride-associated DNA methylation.

BACKGROUND: DNA methylation, an epigenetic modification, can be affected by environmental factors and thus regulate gene expression levels that can lead to alterations of certain phenotypes. Network analysis has been used successfully to discover gene sets that are expressed differently across multiple disease states and suggest possible pathways of disease progression. We applied this framework to compare DNA methylation levels before and after lipid-lowering medication and to identify modules that differ topologically between the two time points, revealing the association between lipid medication and these triglyceride-related methylation sites. METHODS: We performed quality control using beta-mixture quantile normalization on 463,995 cytosine-phosphate-guanine (CpG) sites and deleted problematic sites, resulting in 423,004 probes. We identified 14,850 probes that were nominally associated with triglycerides prior to treatment and performed weighted gene correlation network analysis (WGCNA) to construct pre- and posttreatment methylation networks of these probes. We then applied both WGCNA module preservation and generalized Hamming distance (GHD) to identify modules with topological differences between the pre- and posttreatment. For modules with structural changes between 2 time points, we performed pathway-enrichment analysis to gain further insight into the biological function of the genes from these modules. RESULTS: Six triglyceride-associated modules were identified using pretreatment methylation probes. The same 3 modules were not preserved in posttreatment data using both the module-preservation and the GHD methods. Top-enriched pathways for the 3 differentially methylated modules are sphingolipid signaling pathway, proteoglycans in cancer, and metabolic pathways (p values < 0.005). One module in particular included an enrichment of lipid-related pathways among the top results. CONCLUSIONS: The same 3 modules, which were differentially methylated between pre- and posttreatment, were identified using both WGCNA module-preservation and GHD methods. Pathway analysis revealed that triglyceride-associated modules contain groups of genes that are involved in lipid signaling and metabolism. These 3 modules may provide insight into the effect of fenofibrate on changes in triglyceride levels and these methylation sites.

Correction: Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults.

[This corrects the article DOI: 10.1371/journal.pgen.1006528.].

Genome-wide meta-analysis of 241,258 adults accounting for smoking behaviour identifies novel loci for obesity traits.

Few genome-wide association studies (GWAS) account for environmental exposures, like smoking, potentially impacting the overall trait variance when investigating the genetic contribution to obesity-related traits. Here, we use GWAS data from 51,080 current smokers and 190,178 nonsmokers (87% European descent) to identify loci influencing BMI and central adiposity, measured as waist circumference and waist-to-hip ratio both adjusted for BMI. We identify 23 novel genetic loci, and 9 loci with convincing evidence of gene-smoking interaction (GxSMK) on obesity-related traits. We show consistent direction of effect for all identified loci and significance for 18 novel and for 5 interaction loci in an independent study sample. These loci highlight novel biological functions, including response to oxidative stress, addictive behaviour, and regulatory functions emphasizing the importance of accounting for environment in genetic analyses. Our results suggest that tobacco smoking may alter the genetic susceptibility to overall adiposity and body fat distribution.

Genome-wide physical activity interactions in adiposity - A meta-analysis of 200,452 adults.

Physical activity (PA) may modify the genetic effects that give rise to increased risk of obesity. To identify adiposity loci whose effects are modified by PA, we performed genome-wide interaction meta-analyses of BMI and BMI-adjusted waist circumference and waist-hip ratio from up to 200,452 adults of European (n = 180,423) or other ancestry (n = 20,029). We standardized PA by categorizing it into a dichotomous variable where, on average, 23% of participants were categorized as inactive and 77% as physically active. While we replicate the interaction with PA for the strongest known obesity-risk locus in the FTO gene, of which the effect is attenuated by ~30% in physically active individuals compared to inactive individuals, we do not identify additional loci that are sensitive to PA. In additional genome-wide meta-analyses adjusting for PA and interaction with PA, we identify 11 novel adiposity loci, suggesting that accounting for PA or other environmental factors that contribute to variation in adiposity may facilitate gene discovery.

Visceral Adipose Tissue Is Associated With Bone Microarchitecture in the Framingham Osteoporosis Study.

Obesity has been traditionally considered to protect the skeleton against osteoporosis and fracture. Recently, body fat, specifically visceral adipose tissue (VAT), has been associated with lower bone mineral density (BMD) and increased risk for some types of fractures. We studied VAT and bone microarchitecture in 710 participants (58% women, age 61.3 ± 7.7 years) from the Framingham Offspring cohort to determine whether cortical and trabecular BMD and microarchitecture differ according to the amount of VAT. VAT was measured from CT imaging of the abdomen. Cortical and trabecular BMD and microarchitecture were measured at the distal tibia and radius using high-resolution peripheral quantitative computed tomography (HR-pQCT). We focused on 10 bone parameters: cortical BMD (Ct.BMD), cortical tissue mineral density (Ct.TMD), cortical porosity (Ct.Po), cortical thickness (Ct.Th), cortical bone area fraction (Ct.A/Tt.A), trabecular density (Tb.BMD), trabecular number (Tb.N), trabecular thickness (Tb.Th), total area (Tt.Ar), and failure load (FL) from micro-finite element analysis. We assessed the association between sex-specific quartiles of VAT and BMD, microarchitecture, and strength in all participants and stratified by sex. All analyses were adjusted for age, sex, and in women, menopausal status, then repeated adjusting for body mass index (BMI) or weight. At the radius and tibia, Ct.Th, Ct.A/Tt.A, Tb.BMD, Tb.N, and FL were positively associated with VAT (all p-trend <0.05), but no other associations were statistically significant except for higher levels of cortical porosity with higher VAT in the radius. Most of these associations were only observed in women, and were no longer significant when adjusting for BMI or weight. Higher amounts of VAT are associated with greater BMD and better microstructure of the peripheral skeleton despite some suggestions of significant deleterious changes in cortical measures in the non-weight bearing radius. Associations were no longer significant after adjustment for BMI or weight, suggesting that the effects of VAT may not have a substantial effect on the skeleton independent of BMI or weight. © 2016 American Society for Bone and Mineral Research.